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Abstract 

A recent Faraday wave experiment with two-frequency forcing reports two types of 
'superlattice' patterns that display periodic spatial structures having two separate 
scales [1]. These patterns both arise as secondary states once the primary hexagonal 
pattern becomes unstable. In one of these patterns (so-called 'superlattice- two') the 
original hexagonal symmetry is broken in a subharmonic instability to form a striped 
pattern with a spatial scale increased by a factor of 2\/3 from the original scale of 
the hexagons. In contrast, the time-averaged pattern is periodic on a hexagonal 
lattice with an intermediate spatial scale (\/3 larger than the original scale) and 
apparently has 60° rotation symmetry. We present a symmetry-based approach to 
the analysis of this bifurcation. Taking as our starting point only the observed 
instantaneous symmetry of the superlattice-two pattern presented in [1] and the 
subharmonic nature of the secondary instability, we show (a) that a pattern with 
the same instantaneous symmetries as the superlattice-two pattern can bifurcate 
stably from standing hexagons; (b) that the pattern has a spatio-temporal symmetry 
not reported in [1]; and (c) that this spatio-temporal symmetry accounts for the 
intermediate spatial scale and hexagonal periodicity of the time-averaged pattern, 
but not for the apparent 60° rotation symmetry. The approach is based on general 
techniques that are readily applied to other secondary instabilities of symmetric 
patterns, and does not rely on the primary pattern having small amplitude. 
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1 Introduction 



The classical hydrodynamic problem of parametrically driven surface waves - or Faraday 
waves - concerns the spontaneous generation of standing waves at the free surface of a 
horizontal layer of fluid when subjected to vertical oscillations whose amplitude exceeds 
a critical value. Its usefulness as a tool to study nonlinear pattern-forming dynamics in 
non-equilibrium systems is reflected in the considerable amount of interest shown in the 
subject by experimentalists and theoreticians alike. A review of earlier works, mostly 
conducted with low-viscosity fluids in small vessels and a single forcing frequency, can 
be found in [2]. More recently, Edwards & Fauve [3] have performed experiments in 
the small-depth, high-viscosity and large-aspect ratio regime using a forcing function 
with two commensurate frequency components that modulates gravity periodically. In 
this regime, where it can be shown that the wavenumber of the selected pattern is less 
sensitive to the size and shape of the container, and that long-wavelength modes are 
heavily damped, observations of spatially periodic patterns (stripes, squares, hexagons), 
circular patterns (targets and spirals) and quasi-patterns have been reported [3,4]. A 
survey of more recent results has been carried out by Miiller et al. [5] . Over the past two 
years, a new class of 'superlattice patterns' has been independently observed by KudroUi 
et al. [1] and Arbell & Fineberg [6] in experiments employing two-frequency forcing 
functions, and by Wagner ct al. [7] in experiments using non-Newtonian fluids. These 
superlattices are so termed because of their distinctive feature of having spatial structures 
on two different length scales when viewed at any instant in time [1]. Steady patterns 
that display similar characteristics have also been observed in convection experiments on 
fluids with temperature-dependent viscosity [8] and have been investigated in a model of 
long wave convection [9] and in reaction-diffusion systems near a Turing bifurcation [10]. 

Two types of superlattice patterns have been reported in [1] for different parameter val- 
ues. Both of them, despite their different spatial and spatio-temporal symmetry proper- 
ties, are found to be possible transitions from harmonic standing hexagons as the forcing 
amplitude is increased. (In this context, harmonic indicates an oscillation with the same 
period as that of the external forcing, denoted by T, while subharmonic indicates an 
oscillation with twice that period.) The first of these patterns (called 'superlattice-one' 
by KudroUi et al. [1]) is a harmonic response with triangular symmetry on a small scale 
and hexagonal lattice periodicity on a larger scale. This pattern has been studied by 
Silber & Proctor [11], who showed that it (along with standing hexagons) can arise in a 
bifurcation from the flat, undisturbed state when a hexagonal lattice with spatial peri- 
odicity larger than that dictated by the critical wavelength is considered [12,13]. Silber 
& Proctor [11] also suggested that stabihty might be transfered from standing hexagons 
to superlattice-one through an intermediate branch. 
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Fig. 1. (a) Time-averaged image of the superlattice-two pattern displays a well-deGned hexag- 
onal symmetry on two spatial scales, (b) &: (c) Instantaneous snapshots of the same pattern 
separated by 3/20tii of the external driving period T reveals a time-dependent stripe-like mod- 
ulation. The pattern in (a) is a different realization of the experiment from those in (h) & (c). 
Courtesy of Kudrolli et al. [1], reproduced with permission. 

The second type of superlattice pattern ('superlattice-two' [1]), in contrast to the first, 
arises in a period-doubhng (or subharmonic) instabihty of the standing hexagons. If we 
let m(x, t) measure the deformation of the free fluid surface at time t, it satisfies 

u{x,t + 2T)^u{yi,t), u{x,t + T)y^u{x,t). (1) 

Further, this pattern exhibits a complicated mixture of spatial symmetry and time- 
dependent behaviour. When averaged over two periods of the driving function, its image 
displays hexagonal symmetry with two well-defined spatial scales in the ratio 1 : \/3 
(figure 1(a)). Remarkably, at any instant, a wavy, stripe-like spatial modulation destroys 
the average hexagonal symmetry, resulting in a pattern that appears vastly different from 
its time-averaged image (see figures 1(b) and 1(c)). Arbell & Fineberg (unpublished) have 
also found the superlattice-two pattern for similar experimental parameters. 

The superlattice-two pattern presents a number of theoretical challenges that motivate 
this paper: the disappearance of the stripes from the time-averaged pattern; the reduced 
spatial period in the time-average; and the apparent 60° rotation symmetry of the time- 
average. We present a symmetry-based approach to the study of this pattern by taking 
the view that it arises as a symmetry-breaking instability from the underlying standing 
hexagons in a spatial period-multiplying bifurcation. Our aim is to classify qualitatively 
the range of possible bifurcating solutions and to understand how their symmetry prop- 
erties can be related to the experimental observations described above. We emphasize 
that we are examining instabilities of fully nonlinear states, so our approach differs from 
weakly nonlinear studies of the primary Faraday instability. 
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There are three stages in our approach. First, by using the experimentally observed 
instantaneous spatial symmetry information of the superlattice-two instability and by 
making the assumption that all solutions are periodic in the plane, we can restrict all 
patterns to a suitably chosen spatially periodic lattice. This lattice in turn defines a 
compact symmetry group, which we denote by F*, with the key properties that its 
action leaves standing hexagons invariant and that it has a subgroup, which we denote 
by S'*, that describes the instantaneous spatial symmetry of the bifurcating superlattice 
pattern. Due to the compactness and special structure of F*, we can compute explicitly 
all its irreducible representations. Second, we observe that since is by definition the 
isotropy subgroup of the bifurcating solution under the action of F*, it must have a non- 
trivial fixed-point set. This restriction allows us to identify the one relevant irreducible 
representation of F* that describes the spatial symmetry properties of the marginal 
eigenfunctions at the superlattice bifurcation point. Finally, by considering the action 
of the time-shift symmetry Tj : t ^ t + T on the period-doubling marginal modes, we 
obtain the irreducible representation of the full symmetry group (denoted by F) and 
hence the normal form of the bifurcation problem. We can then invoke the equivariant 
branching lemma [14] to show that there are at least six primary branches of solutions 
bifurcating from standing hexagons. 

With one proviso, the superlattice-two pattern observed by KudroUi et al. [1] can be 
identified as one of these branches, which we show can bifurcate as a stable branch 
from standing hexagons. By applying techniques for studying the averaged symmetry of 
periodic orbits (cf [15]), we show that the time-average of this branch of solutions has the 
hexagonal lattice periodicity observed in the experiment (as in figure 1(a)); this change 
in the spatial length scale on time-averaging is a consequence of the branch of solutions 
possessing a spatio-temporal symmetry. This symmetry was not reported in [1]. The 
proviso mentioned above is that our time-average pattern does not possess 60° rotation 
symmetry - we will return to this discrepancy in the final section. 

A further stability analysis predicts that other patterns, displaying different spatial and 
spatio-temporal symmetry properties, can bifurcate as stable branches of solutions from 
standing hexagons in different regions of parameter space. More generally, our analysis 
indicates that patterns that display superlattice structures can arise in two-dimensional 
spatial period-multiplying bifurcations from an underlying non-trivial solution, and our 
approach could also be used to investigate other superlattice patterns of Arbell & 
Fineberg [6] and Wagner et al. [7]. In particular, it may be possible to analyse some 
of those experimental results in terms of other irreducible representations of the same 
group F*. 

The issue of spatial period-multiplying instabilities is an interesting one that has arisen in 
a variety of experimental and theoretical contexts. Period-multiplying bifurcations in one 

lateral direction have arisen in convection problems [16], magnetoconvection [17], Taylor- 
Couette experiments [18] and in numerical solutions of the Kuramoto-Sivashinsky equa- 
tions [19,20]. Much less is known about spatial period-multiplying bifurcations in two 
directions. There are now several experimental observations of this phenomenon in the 
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Faraday wave problem [1,6,7] as well as in convection experiments [8,16,21] and magne- 
toconvection calculations [22,23]. 

In the next section, we introduce some fundamental definitions and results from equiv- 
ariant bifurcation theory [14] to help us describe how this problem can be cast into a 
theoretical framework. In section 3, we fully describe the symmetry group of the bifurca- 
tion problem that will give rise to the observed symmetry-breaking behaviour. We also 
show that, under suitable phenomenological assumptions, we can identify and hence 
explicitly compute the irreducible representation that is relevant to the action of the 
symmetry group on the observed bifurcating modes. The normal form of the bifurcation 
problem and a stability analysis are presented in section 4. Discussions of our approach 
and a comparison with the experiments follow in section 5. 



2 Group theoretic ideas 

In order to study the superlattice-two pattern as a symmetry-breaking instability from 
standing hexagons, it is necessary to identify all the symmetries that arc initially present. 
Due to the apparent absence of side- wall effects in the observed patterns, we consider the 
mathematical idealisation that all physical fields are defined in a laterally unbounded do- 
main. Standing hexagons are then easily seen to be invariant under the action generated 
by a reflection, a 60° rotation, and two linearly independent translations. The group gen- 
erated by these symmetry actions is isomorphic to -i-Dg, which is non-compact. (Here 
Z denotes the group of integers under addition, and Dg is the twelve-element symmetry 
group of a regular hexagon.) Consequently a bifurcation problem that is equivariant un- 
der the action of this group can have an infinite number of modes related by symmetry 
becoming marginally stable simultaneously. This difficulty can be resolved if we restrict 
possible solutions to doubly-periodic functions defined on a suitably chosen lattice, an 
assumption justified by the distinct spatial periodicity of the observed patterns. A suit- 
able lattice, which can be viewed as a finite cell with periodic boundary conditions, is 
one that captures the spatial periodicity of both the bifurcating modes as well as the 
standing hexagons. With respect to such a periodic cell, the symmetries that leave stand- 
ing hexagons invariant now form a finite, and hence compact group that can be studied 
via representation theory. Our task therefore, is to make use of the available symmetry 
information taken from experimental observations to choose a lattice on which we can 
define a suitable spatial symmetry group with the properties outlined in section 1. 
The idea of suitability can be made precise after we have introduced some basic group 
theoretic results. 

Since we are considering bifurcations from a time-periodic solution, we formulate the 
bifurcation problem of the superlattice-two pattern, a period-doubling instability, by ex- 
panding about standing hexagons using a stroboscopic map Q in the manner described 
by Crawford & Knobloch [24] and Silber & Proctor [11]. Specifically, we are assuming 
that standing hexagons, a fixed point of the map Q, lose stability to subharmonic waves 
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with period 2T as a bifurcation parameter /i is varied past zero. This imphes that the 
hnearised map DQ evaluated at the fixed point has a real eigenvalue passing through the 
value —1 as the bifurcating waves become unstable. With the assumption that all fields 
are defined in a periodic cell such that symmetries of the standing hexagons are described 
by a compact group F*, the linearised map DQ has a finite number (p) of marginal eigen- 
functions associated with the eigenvalue — 1 as ;U crosses zero. We denote the amplitudes 
of these p marginal modes at time t = qT, g G Z by = [zi{qT), . . . ,Zp{qT)] G MP. 
In addition, the pattern has two neutrally stable modes (eigenvalues equal to 1) asso- 
ciated with translations of the standing hexagons (see [18,25]); the amplitudes of these 
two modes, which correspond to the translation of the pattern in the plane, are de- 
noted by dq. Close to the onset of the period-doubling instability, Q can be reduced to a 
finite-dimensional map g defined on the centre manifold spanned by these amplitudes: 

z,+i = g(z,;//), grR^'xR^R^ (2) 

coupled with a map h : R^ x R — > R^ describing how the perturbation drives translations 
of the pattern: 

dq+i = d5 + h(z5;/x). (3) 
The map g is forced by symmetry to be F^'-equivariant: 

7g(zg; = 8(72?; for all 7 e F^ (4) 

while the map h obeys 

N^h {zg-, 11) = h {jzg] 11) for aU 7 e F*, (5) 

where is the 2x2 matrix that represents how the symmetry 7 acts on a horizontal 
displacement vector [25]. In terms of the marginal modes, standing hexagons correspond 
to the trivial state, z = 0. When considered as a symmetry-breaking bifurcation from 
the underlying standing hexagons, the superlattice pattern corresponds to a non-trivial, 
period- two solution to the map g, denoted by z*, whose instantaneous spatial symmetry 
is specified by its isotropy subgroup E^* : 

^{aer-. azl = z;} C F^ (6) 
In fact, this solution must lie in the fixed-point subspace of E^. : 

Fix(E^*) ^{zeW -.az^z, for all (7 e E^ J , (7) 
which is a linear subspace of R^ and invariant under g [14] . 

Since standing hexagons do not possess spatio-temporal symmetries [26], F*-equivariance 
is sufficient to determine the normal form of the map g in the case of (temporal) period- 
preserving bifurcations. However, for period-doubling bifurcations there is an extra sym- 
metry pertaining to the normal form, related to the time-shift action Tt : t ^ t + T 
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on the bifurcating modes. In this case, the normal form of the map g is (F* x Z2)- 
equivariant [27,28] . Once the full symmetry group of the normal form of g is determined, 
we can apply the equivariant branching lemma, which, with suitable interpretation, 
states that if certain non-degeneracy conditions are satisfied, there is a unique branch of 
bifurcating solutions for each isotropy subgroup of F = x Z2 with a one-dimensional 
fixed-point subspace. So instead of solving for solutions of the nonlinear vector field g, we 
can simply look for isotropy subgroups of F with this property. To apply the equivariant 
branching lemma, we need to know explicitly how symmetry acts on all the marginal 
modes, but experimental observations only provide information about the instantaneous 
spatial symmetry of one of these modes. We cannot infer directly from the observations 
the total number of marginal modes that are related by symmetry at the bifurcation 
point, nor the set of matrices that represent the action of the symmetry group F on the 
marginal modes and the map g. However, this difficulty can be resolved if we make the 
(generic) assumption that the bifurcation is associated with an irreducible representation 
of the group F, and this is where the need to invoke representation theory arises. 

In order to introduce the key properties of irreducible representations (irreps) and de- 
scribe how they can be computed for a finite group F, we recall the following defini- 
tions [29]. 

(i) A representation of the group F is a homomorphism ip that maps F into a set of 
invertible n x n matrices Mp acting on M" or C", in other words 

^ (7) ^M^, 7 e F, e Mr 

such that (7172) = 4'{li)ip{l2) for all 71, 72 € F. The integer n is the dimension of 
the representation. 

(ii) Two n-dimensional representations Mr and Nr of F are called equivalent if there is 
an invertible n x n matrix Q such that for each 7 G F, 

= Q-^M^Q, e Mr, G Nr. 

(iii) A conjugacy class of F is a subset C of F such that 7~^C7 e C, for all c e C and 7 e 
F. 

(iv) The character of an element 7 G F in a representation Mr is defined to be the trace 
of the matrix M^, and we denote this value by Xm^- 

(v) A representation of F on M" (C" ) is said to be irreducible if it does not leave invariant 
any proper subspace of R" (C"). 

Simplistically we can consider a representation as a set of n x n nonsingular matrices 
that specifies the action of F on the vector space R" or C" and at the same time preserves 
the group structure. It is possible to show that every representation of a finite group 
is equivalent to a unitary representation - one in which all matrices are unitary [29]. 
A simple result of definition (iv) is that the character of the identity element in a 
representation is always equal to the dimension n of that representation, and definitions 
(i), (iii) and (iv) imply that elements in the same conjugacy class have the same character. 



7 



The characters of the irreps of F obey a set of rules inherited from the orthogonaUty 
theorem governing the underlying irreps [30] and for simple groups such as Z2, Zq and 
Dg, the character tables can easily be constructed by appealing to those rules. The 
orthogonality theorem also implies that the number of irreps of a group is equal to the 
number of conjugacy classes. For finite groups with a semi-direct product structure of 
the form F = ^ -j- B such that ^ is a normal (or invariant) subgroup of F (that is, 
7~^a7 G A for every a & A and 7 G F), the characters of the irreps of A and B form 
the building blocks in determining all the characters and constructing unitary irreps of 
the group F via a special algorithm [29] . 

In summary, analysis of the superlattice pattern using these group theoretic tools de- 
pends on our being able to find a spatial lattice or periodic cell on which standing 
hexagons and the marginal modes exhibiting the observed symmetries fit. The arrange- 
ment of the standing hexagons in the periodic cell then gives us a suitable symmetry 
group F*, which has a subgroup E^. , defined in (6), whose elements are determined from 
experimental observations. Once we have calculated all the characters of F^, the restric- 
tion provided by the requirement that S^. be the isotropy subgroup of the observed 
pattern enables us to isolate the one irrep that describes the action of F* on all the 
marginal modes related to the observed bifurcating mode by symmetry. Indeterminacy 
in the choice of irreps can be avoided if we choose a unit cell that captures exactly one 
spatial period of the observed pattern. The details of this procedure are the subject of 
the next section. 




Fig. 2. (a) Schematic representation of standing hexagons, where ei and 62 denote the vectors 
of translations defined in (10). (b) Diagram depicting the instantaneous spatial symmetry 
of the superlattice-two pattern (see figure 1(c)), whose spatial periodicity can be captured 
by hexagonal cells mapped to one another by the translations ei and 62 as shown in (c). 
The superlattice-two pattern is left unchanged by r|, Kx and rfp^. The small circles, dotted 
lines and shading in (b) and (c) serve to identify equivalent hexagons related by translations, 
rotations and reflections of the pattern. 
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3 Finding the symmetry group F and the irrep for the superlattice-two 
bifurcation problem 



A closer examination of images obtained from the experiment reveals that it is possible 
to impose a hexagonal lattice on the observed patterns, whose instantaneous spatial 
symmetries are depicted in figure 2. The choice of lattice is not unique as it can be shown 
that there are many possible candidates (for example a \/3 : 1 rectangular lattice), but 
a hexagonal lattice is a natural choice due to the symmetry of the standing hexagons. 
Let us denote the two generating vectors of the hexagonal lattice £ by ei, e2 G such 
that 

|ei| = |e2|=c, (8) 

where c is a scaling factor (figure 2(c)). Functions in the plane that are doubly-periodic 
with respect to £ satisfy 

u{x,t) ^u{x + \,t), x= (x,y) e K^ 1 e £, (9) 

where the lattice is defined as 

C = jniei + n2e2 : (ni, 712) G Z^j . 



First let us consider the spatial symmetries of the standing hexagons shown schematically 
in figure 2(a). They are invariant under the action of Dg as well as two translations, which 
we define as follows: 

Ti:x^x + ei, r2:x— >x + e2, (10) 

and let |ei| = |e2| = cq be the observed size of the periodic cell in which the basic 
standing hexagons fit. Our aim is to pick a value for the scaling factor c in (8) in terms 
of Cq. As indicated at the end of section 2, a suitable choice of the value c is one that 
gives a hexagonal cell whose size captures precisely one spatial period of the bifurcating 
modes, as shown schematically in figure 2(c). In fact, the observed ratio of the two 
lengths |ei| and |ei| is c/cq = 2-\/3, and for this value of c, the symmetry group F* 
of the standing hexagons includes non-trivial translations generated by Ti and T2. The 
structure of the group {ti,T2) can be determined if we express each of the translations 
in (10) in terms of Ci and 62. We can then look for the lowest powers ni,n2, n^, n4 G Z+ 
such that Ti\ r^^, rf'V^'' map the lattice C to itself, and thus determine the order of 
the group (ti,T2). 

Guided by the experimental observations, wc choose ei = Cq (^^, |)i ^2 = cq (0, 1) such 

that ei = 4ei - 2e2 = co{2V3,0), 62 = 2ei + 2e2 = Co(V3,3) (see figure 2). The 
translations can now be written as 

n : x-F -ei -e2, T2 : x ^ x - -ei -62, (11) 
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and we can easily show that they satisfy = r| = r| = identity, as vectors of the 
form X + TTT-iei + 771262 for any integers mi and m2 lie in £ and arc therefore identified. 
Since Ti and T2 commute we can also see that every element generated by Ti and T2 can 
be written as t"t2 or rf for n = 0, . . . , 5. In total there are twelve different translations, 
forming a group that is isomorphic to Zq x Z2. The order of this group being twelve 
corresponds to the fact that each of the large hexagonal cells in figure 2(c) contains 
exactly twelve of the smaller hexagons. 
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Table 1 

Character table of the group of spatial symmetries constructed via the algorithm taken 
from [29]. A representative element is shown for each conjugacy class, and the number of 

elements in the class is given in brackets. Classes marked by * contain elements (specified at 
the top of the table) of the eight-element group E^* = ( tI , TiP^)- 



So in terms of the lattice £, the full spatial symmetry of the standing hexagons is given 
by the group = {Zq x Z2) + De, where De is generated by a refiection k.^ and a 
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60° rotation p and its standard action on is given by 

i^x-{x,y)^{-x,y), p : (x, ^ ^ (x - VSy, VSx + , (12) 
and Zg X Z2, an invariant subgroup of F*, is generated by the two translations 

Ti : (a;,^/) ^ + ^co,?/ + ^coj , T2 : {x,y) ^ {x,y + cq) . (13) 

The group has the semi-direct product structure mentioned in section 2. As a result 
we can apply the algorithm taken from [29] to calculate all its characters and irreps, and 
we present the characters of its irreps in table 1. 

Any elements that have the same character as the identity in a unitary irrep of must 
also act like the identity [31]. Using this simple idea and the information taken from 
experimental observations about the spatial symmetries of the unstable mode, we can 
single out the irrep of that describes the instantaneous symmetry-breaking behaviour. 
Careful study of snapshots of the supcrlattice-two pattern shows that it is invariant under 
the action of r|, Kx and rfp'^ (sec figures 1(b) and 1(c), where the pattern is shown at 
a slightly tilted angle, and figure 2(b)). The group generated by these elements is by 
definition the isotropy subgroup of the bifurcating mode under the action of F* (cf (6)), 
therefore 

^(r|, riV)cF^ (14) 

We can now go through the hst of characters of F* given in table 1 and determine which 
irrep satisfies the criteria of permitting E^* defined in (14) to be an isotropy subgroup 
of the bifurcating solution. First, any irreps that satisfy 

Xm^, = XMia for some 7^ e F" and 7^ ^ E^. (15) 

must be rejected, because in these cases the isotropy subgroup of z* must contain spatial 
symmetry elements apart from those that are observed. This eliminates representations 
1-12. In representation 13, the class containing r| is represented by the identity, but this 
class also contains rf and T1T2, which are not in E^., so eliminating this irrep and leaving 
only 14 and 15. Then we can use the trace formula [14] to calculate the dimension of 
the fixed-point subspace of E^. : 

dim Fix (E) = ^ ^ xm,, 



which gives for representation 14 and 1 for representation 15. Clearly we require 
dim Fix (^E^.^ 7^ 0, since Fix(^E^*^ is non-trivial. Thus the six-dimensional irrep M^s 
is the only one in which E^* satisfies the conditions of being an isotropy subgroup of the 
observed mode. 
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In addition to being equivariant under the action of spatial symmetries as specified by 
this irrep, the normal form of the period-doubling bifurcation problem has an extra 
symmetry corresponding to a translation in time by one period of the external forcing: 

Tt-.t^t + T. (16) 

This element can be viewed as a spatio-temporal symmetry with a trivial spatial action, 
and it acts independently from elements in with respect to the standing hexagons. So 
the full symmetry group F of the normal form for the superlattice bifurcation problem is 
a direct product between F* and the group (rt), which, as can be seen from (1) and (16), 
is isomorphic to Z2, hence F = F* x Z2 as we pointed out in section 2. We can write each 
element 7 e F as 

7 = (7., <7t) , 7. e ^^ at e (n) , (17) 

such that for 71 = (7^1, cxtj, 72 = (7^2, cxtj, 7172 = {'Jsilsi, (^tif^ti)- Because of the direct 
product structure of F and the period-doubling nature of the bifurcating solution, 
must act like —1 on the amplitudes of the marginal modes. Therefore the irrep of F that 
specifies the action of spatial and spatio-temporal symmetry elements on the marginal 
modes and the normal form of g can be constructed from the set of matrices My^ e Mps 
as follows: 

^ M-y^ if at — identity 
^ -M^^ if at = Tt 

for all 7 = (7s, at) G F. This irrep, which we denote by Mr, is of the same dimension 
as Mps , which implies that we have a six-dimensional centre manifold at the bifurcation 
point. So all bifurcating solutions can be written as m(x, t) = uo{x, t) + C(x, t) such that 

C(x, qT) = AJi(x) + Sj2(x) + C,/3(x) + c.c. + h.o.t., q G Z (18) 

where uo{x,t) represents standing hexagons, c.c. denotes complex conjugate, h.o.t. de- 
notes the higher-order terms, and Ag, Bg, Cg G C are the small amplitudes of /i, /2 
and /a, the three complex marginal eigenfunctions that form a basis for the neutral 
eigenspace (excluding the two zero eigenvalues corresponding to translating the under- 
lying pattern). Note that by including the higher-order terms, ^ represents the nonlinear 
perturbation from the standing hexagons. 

Applying the method described in [29], we can construct all the 6x6 matrices 
that specify the action of F on (or C'^) for the irrep Mp. Rather than describe this 
procedure, we find it convenient to specify the group action by choosing a small number 
of Fourier modes to represent the marginal eigenfunctions, and working out how the 
amplitudes of these modes Ag, Bg and Cg transform under the generating elements of F. 
Since representations are defined only up to a similarity transformation, the choice of 
Fourier modes we make will not matter, as long as we are careful not to introduce any 
accidental symmetries (which would become apparent on checking the characters). 
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Any function w(x, t) defined on the lattice C can be written as a double Fourier series 
of the form 



^) = E E (t) e'-^i^^^^+^^^^y-, (19) 

where ki and k2 are the generating wavcvcctors of the dual lattice C* related to ei and 
62 by kj ■ = 6ij such that (9) holds. Our choice of the vectors ei and 62 requires 
ki = /c and k2 = A;(0, 1), where k ^ 

We use the observed instantaneous symmetry of the pattern (see figure 2(c)) to select a 
representative function from the full set of Fourier modes, starting with a single Fourier 
mode (>'^'^^ih^-i+h^2)-^ (and its complex conjugate) for some choice of integers ji and j2- 
If the pattern is to be invariant under r|, ji must be even, so set ji = 2m, and, for later 
convenience, set j2 = m + n, where m and n are integers. With this choice, the Fourier 
mode is e^'^''^(^"*^+"^). The pattern is also invariant under rfp^. Now replaces the 
chosen mode by its complex conjugate, and rf multiplies the mode by a complex number 
with unit modulus. Since rf is of order two but not equal to the identity, it must act 
by multiplying the mode by —1. This forces m + n to be odd (and so for the observed 
pattern, the amplitude of the Fourier mode must be pure imaginary). The translation ri 
must act with order 6 (otherwise the pattern would be invariant under a lesser translation 
in that direction), so 3m + n = 1 mod 6 or 3m + n = 5 mod 6. The second of these is 
essentially the complex conjugate of the first, so we choose 3m + n = 1 mod 6; (m, n) 
could be (0,1), (2,1) or (1,4), for example. Finally, the reflection generates a new 
function ^'^'^^H-^'^^+i'iy) ^ so the superlattice-two pattern can be exemplified by a mode 
of the form /i = e^nikiVsmx+ny) ^2nik{-V3mx+ny) _ gj^^y degree rotations of this function 
generate /2 and /a, so we have: 

_ g27riKi-x _|_ g27riK2-x _ g27riK3-x _|_ ^2mK4-x _ ^2TviKs-x _|_ ^2TviKe-x (20) 

(the true eigenfunctions will be made up of linear combinations of such functions) , where 

Ki = A;(v^m, n), K.2 = k{—^/3m,n), 

Ks = I (V3(m + n),(-3m + n)), K4 = | (V3(-m + n), (3m + n)), 

K5 = I (V3(m + n), (3m - n)). Kg = | (v^(-m + n), -(3m + n)). 



These wavevectors have the same wavenumber K{m, n) = k^/3m? + n?, with m and n 
satisfying 3m + n = 1 mod 6. With this choice of basis functions, the relevant irrep 
of r can be specified by the action of the generating elements of F on the amplitudes 

(Aq, Bq, Cq): 
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K-x '■ i^q, Bq, Cq) — > {Ag, Cq, Bq), (21) 

p:(Aq,Bq,Cq)^{Bq,Cq,Aq), (22) 

Tl : {Aq,Bq,Cq)^{efAq,e-^Bq,efCq), (23) 

T2 : {Aq, Bq, Cq) ^ Aq, 6^ Bq, C'^^Cq), (24) 

Tt : (A„ S„ Cq) ^ (-A„ -C,). (25) 



We include the subharmonic action of Tt here for completeness. The same representation 
could be constructed using the method described in [29]. 

4 Normal form of the bifurcation problem 

We now have sufficient information to invoke the equivariant branching lemma [14] 
and describe the different patterns that must be formed in the instability that created 
the superlattice-two from standing hexagons. Before doing this, we will compute the 
normal form for the bifurcation since we need it to work out the stability of the various 
patterns. The irrep (21-25) we identified in section 3 implies that the reduced map g 
introduced in (2) is six-dimensional, and we let Zq — {Aq, Bq, Cq), g G Z, Aq, Bq, Cq e C. 
As indicated earlier, the action of Tt defined in (16) is due to the subharmonic nature 
of the bifurcating modes with respect to the overall driving period T given in (1). If 
each iteration in Zq corresponds to advancing in time by T, then = —'^q+i = Zg. 
Consequently, g(zg) = z^+i = -Zq+2 = -g(zq+i) = -g(-Zq) [32]. So the map g will be 
an odd function of the amplitudes Ag, Bg and Cg, as well as being F-equivariant. This 
information enables us to write down the form of g including up to fifth order terms: 

Aq+i = - (1 + /x) Ag + ailA/A, + «2 {\Bgf + \Cgf) Aq + Al 

+ + \CX) Ag + /3s\Agf {\Bqf + \Cqf) Ag + /3,\ B g^ \C g^ Ag 

+ ^.B'^C'gAq + uA^, 

Bq+i ^-{l + n)Bq + ai\BqfBq + (}Aq\^ + |C/) Bq + 

+ (52 (1^1/ + |C/) Bg + (3^\Bgf [\Agf ^ \C g^) Bg ^ (5 ^\ Ag^ \C g^ B g 

+ P,AlC'gBg + uBl, 

Cg+, = -{l+^)Cg + a,\CgfCg + 0^2 (| + l ) Cg + Pl\Cg\'Cg 

+ P2 {\Ag\^ + W) Cg + /3s\Cgf {\Agf + \ B g^) Cg + /3,\ Aq^ \ B q^ C g 

+ P,AlBpg + uC'g, 

where all coefficients are forced by symmetry to be real. 

Apart from the u terms, the equations above are equivalent to the T^-i-D6xZ2-equivariant 
amplitude equations (truncated to the same order) that arise in the context of Boussinesq 
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(26) 
(27) 
(28) 



convection on a hexagonal lattice [14,33], once they arc rc-intcrprctcd as amplitude 
equations rather than a map. The v terms have the effect of breaking the full (two- 
torus) symmetry group of translations in a periodic domain to the discrete translations 
allowed by the underlying pattern. A natural question to ask is why we needed to work 
out the details of the representation before writing down these amplitude equations. 
The main reason is that we did not know in advance how many linearly independent 
marginal eigenfunctions are involved in the instability. Even if we had assumed that 
there were six, it has turned out that there are two six- dimensional irreps, only one 
of which is involved in the bifurcation. The other six-dimensional irrcp is generated by 
taking /i = ^2mk{V3mx+ny) _ g2mk{-V3mx+ny) g^j^^ (following a similar analysis) results 

in the same amplitude equation. Without realising this, one might conclude incorrectly 
that patterns that are odd under reflection might also be found in this instability. All 
the other irreps in table 1 have dimension less than six (that is, there are fewer than six 
independent marginal eigenfunctions), so the order of the relevant normal forms would 
be correspondingly less. 

We also use (5) to write down the dynamics of the position of the underlying standing 
hexagons, truncated to quartic order: 



dq+i ^dq + ^lm 



(29) 



where ^ is a constant. 



representative solution branch isotropy subgroup averaged symmetry 


I 


1. ^gM, B = C = (r|, K^, p^, tI) 

2. Aem,B = C = Q (t|, Kx, nV'^ rf) 


(''"i) ^ii '^xt p^) 


II 


3. A = B = C D6 = {nx, p) 

4. A = -B = C {kx, p) 




III 


5. A = 0, B = C eR {kx, p^, rf ) 

6. A = 0, B = -C GiR {kx, tIp^, r|) 


(r^, Hx') P^) 



Table 2 

Primary solution branches of the normal form (26-28) and their isotropy subgroups, grouped 
into three types, I-III. Using the analysis presented in section 4.2 we can show that the 
two branches within each type of solution share the same time-averaged spatial symmetries. 
Branch 2 of type I corresponds to the superlattice-two pattern. 



We can show that there are six isotropy subgroups whose fixed-point subspaces are 
one-dimensional, so the equivariant branching lemma tells us that there are at least six 
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primary bifurcating branches of solutions from standing hexagons, and we summarise 
these solutions in table 2. Elements accented by a tilde represent spatio-temporal sym- 
metries, which, using the notation introduced in (17), can be written as ff = (rf , Tj), and 
similarly for f| and p. The super lattice-two pattern corresponds to branch 2 of type I. 



(a) standing hexagons (b) branch 3 (c) branch 4, t = (d) branch 4, t = T 




(i) branch 5, t = (j) branch 5, t = T (k) branch 6, t = (1) branch 6, t = T 




Fig. 3. Instantaneous planforms of the different solution branches summarised in table 2 and 
illustrated here in frames (b)-(l) as small-amplitude perturbations to standing hexagons. Solid 
squares represent lattice points of C. (a) Standing hexagons, which have the full F symmetry, 
(b) Solution branch 3 with Dg symmetry, referred to in [13] as 'supcrhexagons'. The periodic 
hexagonal boxes are delineated by light borders surrounding each cell, (c) &: (d) Solution branch 
4 att = and T showing D3 symmetry as well as the spatio-temporal symmetry p. (e) &: (f) The 
superlattice-two pattern corresponds to solution branch 2 as they share the same instantaneous 
spatial symmetries. This pattern is shown here at t = and T with spatio-temporal symmetry 
rf evident, (g) & (h) Similar spatio-temporal symmetry is displayed by solution branch 1. 
(i)-(l) Branches 5 &: 6 have very similar symmetry properties: both are invariant under the 
action of rf and instantaneously they differ only by a shift of the reflection symmetry Ky . 
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For the choice of wave integer pair {m,n) = (2,1), the instantaneous planforms of 
these six solution branches are illustrated schematically in figure 3. We can compare 
figures 3(e) and 3(f) with 1(b) and 1(c) and notice that the appearance of stripes at 
regular intervals in the grey-scale plots of solution branch 2 closely resembles the essential 
features of the experimentally observed superlattice-two pattern. 

None of these primary branches leads to a net drift of the underlying hexagonal pattern; 
this can be seen in two ways: first, because the rate of drift (from (29), truncated to 
quartic order) is zero on all six primary branches; second (and more convincing) since 
is in the symmetry group of all the time averaged patterns (see below). In other words, 
the patterns are all pinned by the 180° rotation symmetry on average. 



4-1 Stability results 

We summarise in table 3 the branching equations and the Floquet multipliers of the 
period 2T patterns, for each of the six primary solutions guaranteed to exist by the 
equivariant branching lemma. Floquet multipliers greater than one in magnitude indicate 
instability. We group the six branches into three types and denote them by I, II and III 
as shown in tables 2 and 3. It is evident that branches within each type are degenerate up 
to third-order terms, thus necessitating the inclusion of quintic terms. In particular, only 
one solution branch within each of types I and II can be stable depending on the signs 
of and -|- 1^, and both branches in type III are always unstable. Only one branch can 
bifurcate stably, and all branches must be supercritical for one of them to be stable. One 
of the requirements for the observed superlattice-two pattern (i.e., branch 2 of type I) 
to be stable is that the quintic coefficient z/ > 0. If we also assume the non- degeneracy 
conditions ai ^ 0, ai + 2a-2 7^ 0, ai ± 0:2 7^ 0, z/ 7^ and /Ss + z/ 7^ 0, close to the 
bifurcation point the relative stability of branches of the three types is illustrated by the 
bifurcation diagrams shown in figure 4. 

The experimental results [1] suggest that the bifurcation may have subcritical branches 
as there is a parameter regime in which standing hexagons and the superlattice-two 
pattern may coexist. On the other hand, the experimentalists report no hysteresis be- 
tween standing hexagons and superlattice-two, while they do report hysteresis between 
hexagons and other patterns at other parameter values, so it is not clear whether or not 
there is a direct bifurcation from standing hexagons to the superlattice-two pattern in 
the experiments. With our parameters, we require ai > 0, «2 > and z/ > for the 
superlattice-two pattern (branch 2 of type I) to bifurcate stably, but the branch could 
also be stable in the region cti < 0, q;2 > ai and > if there were a saddle- node 
bifurcation on branch I. 
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Primary solutions and branching equations 


Floquet multipliers (multiplicity) 


I 


1. Aq = Ri, Bq = Cq = 0, 

= + aiRl + (A + u) R\ 


1 - AaiRl, 1 - 2 (0:2 - ai) Rj (4 times), 
1 + 12uRf 


2. Aq = iR2, Bq = Cq = 0, 

= -/X + aiRl + - v) Rl 


1 - AaiR^, 1 - 2 (q;2 - ai) (4 times), 
1 - 12uR^ 


II 


3. Aq = Bq = Cq = R3, 

= -/X + (ai + 2^2) Rl 

+ 2/32 + 2/33 + /34 + /35 + i^) RI 


l-A{ai + 2a2)Rl 

1 + 4 (q;2 - ai) i?| (2 times), 

1 + 12uRl (2 times), 1 + 12 (/35 + i/) R^ 


4. = -Bq = Cq= iR4, 

= -/X + (ai + 2q;2) Rj 

+ (/3i + 2/32 + 2/33 + /34 - /35 - I/) 


1 -4(ai + 2a2)i?i 

1 + 4 (q;2 - ai) (2 times), 

1 - UiyR'l (2 times), 1 - 12 (/35 + i^) Rf 


III 


5. Aq = 0, Bq = Cq = R5, 

= -/X + (ai + 02) Rl 

+ {/3i + p2 + + 1^) Rt 


1 - 4 (ai + a2) iili, 1 + 4 (^2 - ai) itlf, 
1 - 2 (a2 - ai) (2 times), 
1 + 12i/i?| (2 times) 


6. = 0, = -Cq = iRe, 

= -fi + {ai + ^2) Rj 
+ + p2 + - i^) Ri 


1 - 4 (ai + 02) i^G' 1 + 4 ("2 - ai) i?6' 

1 - 2 (02 - cti) (2 times), 
1 - UuR^ (2 times) 



Table 3 

Branching equations and Floquet multipliers for the six primary period-two solutions of the 
normal form (26-28) listed in table 2. Aq, Bq and Cq are complex amplitudes of the marginal 
modes defined in (18). Only leading order terms in R^ are shown, and the multiplicities of the 
Floquet multipliers (computed from the second iterate of the map) are indicated. 



4-2 Time-averaged behaviour 

We can study the symmetry properties of the time- averaged image of the observed solu- 
tion by integrating over a full period of the newly created periodic orbit (cf [15,25,34]). 
Specifically, we let UQ(x,t) be the standing hexagons solution and C(x, t) the nonlinear 
perturbation to this solution such that u{'x, t) = 'Uo(x, t) +C(x, t) represents the observed 
pattern. We know that uo(x,t) is F-invariant, i.e., 7Mo(x, t) = mo(x, t) for all 7 G F, and 
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Fig. 4. Bifurcation diagrams for the T-equivariant normal form (26-28). The sign of the cubic 
coefRcient ai determines whether solution type I bifurcates sub- or supercritically. If we assume 
that the quintic coefRcient u > 0, then the superlattice-two pattern (branch 2 of type I) can oc- 
cur as a stable branch in the region of the (ai, 0:2) space given by {(ai, 0:2) : ai > 0, 02 > ai}. 

have also found that the spatial and spatio-temporal symmetry of C(x, t) is given by 
its isotropy subgroup (t|, Kx, Tfp^, ff) = C F. Let 7^ and 7^ denote respectively 
the purely spatial symmetry elements and the spatial part of spatio-temporal symmetry 
elements in that act on C(x, as follows: 



The time-averaged value of u(x, t) can be obtained by integrating over the full period of 
the bifurcating solution: 



where we have used the fact that Mo(x, t) = Mo(x, t + T). Clearly, u shares the same 
spatial symmetry with ( because both Uq and the individual entries of the integrand in 
(30) are invariant under the action of 7^ G E^. It is also invariant under 7^ because the 
integrand in (30) as a whole is invariant under 7^: 



7,C(x, t) = C(x, t), 7iC(x, t) = C(x, t + T). 




1 



•2T 



(30) 
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_ _ 1 

7tM(x) = 7t'Uo(x) + — / 7tC(x, t) + 7iC(x, t + T)dt 
21 Jo 

= ^o(x) + 2^ C(x, t + T) + C(x, dt 
= fZ(x). 

This result in fact follows readily from more general results on the symmetries of chaotic 
attractors [15,34]. 

In the case of the observed sTiperlattice-two pattern with isotropy subgroup S^, the 
spatial component of the spatio-temporal symmetry clement, namely rf, will show up 
alongside r|, and rfp"^ in the time- averaged image to generate an augmented spatial 
symmetry group = {Tf,T2,nx,p^)- This prediction is in agreement with experimental 
results (figure 1(a)) and can be understood in the following way. The action of the 
translations Ti and T2 on ■u(x) is of order three since rfu{'x.) = r|-u(x) = 'u(x), whereas 
the order of the same action on m(x) is six. So the averaged pattern will appear to be 
periodic on a lattice £av spanned by basis vectors such that |eav| = ||ej| = |c. We 
have shown in section 3 that cq = it follows that |eav| = \/3|ei|. Therefore the 
ratio of spatial period of the averaged pattern to that of the basic standing hexagons is 
1 : which is consistent with the observation reported by [1] as shown in figure 1(a). 
Using the same reasoning and information from the isotropy subgroups of the primary 
solutions given in table 2, we therefore predict both branches in each type of solutions 
to have the same time-averaged symmetries. 



5 Discussion 

Starting from the observed instantaneous symmetry of the superlattice-two pattern re- 
ported in [1] , we have been able to show (a) that a pattern with the same instantaneous 
spatial symmetry as the superlattice-two pattern can bifurcate stably from standing 
hexagons in a spatial period-multiplying instability; (b) that the pattern has the spatio- 
temporal symmetry (not reported in [1]) of advancing one driving period in time com- 
bined with a translation by three units in space (figure 3(e) and 3(f)); and (c) that this 
spatio-temporal symmetry accounts for the intermediate spatial scale and periodicity 
on a hexagonal lattice of the time-averaged pattern (figure 1(a)). We should emphasise 
that the intermediate spatial periodicity of the time-averaged pattern is not the spatial 
periodicity of the larger hexagonal lattice that we have assumed. 

Arbell & Fineberg (unpublished) have found the superlattice-two state in their exper- 
iments and have confirmed that it does have the spatio-temporal symmetry that we 
predict. Our results also suggest that 60° rotations are not in fact symmetries of the 
time-averaged pattern, but should be weakly broken. The breaking of 60° rotational 
symmetry, if it is present, is evidently a small effect since the hexagons in figure 1(a) 
do appear to be invariant under 60° rotations [1] (this has been confirmed by GoUub, 
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private communication). For other parameter values, the symmetry breaking effect may 
be more pronounced: Fineberg (private communication) reports that his experimental 
time-averaged pattern is not invariant under 60° rotations. Clearly this would be an 
interesting issue to investigate in more detail, but the measurements are delicate and 
are liable to be prone to systematic errors or imperfections, so confirming our prediction 
could be difficult. 

The spatio-temporal symmetry of superlattice-two arises because the instability of stand- 
ing hexagons is subharmonic. Other patterns, with different combinations of spatial and 
spatio-temporal symmetries, are possible stable branches in the same bifurcation prob- 
lem. Not all branches of solutions have spatio-temporal symmetries, and some of the 
patterns share the same time-averaged symmetries even though they have different in- 
stantaneous planforms. The method we have presented is based entirely on symmetry 
arguments and is able to deal with instabilities of a fully nonlinear time-periodic solution. 

Spatial period-multiplying instabilities have arisen in a variety of contexts, in both 
one [17-20] and two lateral directions [1,6-8,16,21-23]. Most of these situations involved 
relatively simple groups; part of the difficulty and interest here has been the size of the 
symmetry group, enlarged because of the number of translations broken by the new pat- 
tern. Only one of the 15 representations is involved in the superlattice-two bifurcation; 
other representations may be relevant to other experiments (particularly [6,7]) in which 
standing hexagons lose stability to patterns that fit into the larger hexagonal cells we 
have used here. 

As can be seen in section 3, a heuristic step in our method involves the choice of a 
suitable periodic cell that accommodates the observed patterns and whose size coin- 
cides with exactly one spatial period of the bifurcating modes. The arrangement of the 
underlying basic state in this cell then defines a spatial symmetry group F* of the bifur- 
cation problem and the instantaneous symmetries of the superlattice instability form its 
isotropy subgroup E*. If a larger hexagonal periodic cell that captures more than one 
spatial period of the bifurcating modes had been chosen, the translations Ti and T2 given 
in (10) that leave standing hexagons invariant would have had higher order, resulting in 
a larger spatial symmetry group. In this case there would have been more than one irrep 
of F'* in which satisfied the conditions of being an isotropy subgroup. By choosing the 
smallest possible periodic cell, we have found that such indeterminacy can be avoided. 

The method we have described in this paper for analysing certain types of symmetry- 
breaking instabilities bifurcating from a non-trivial basic state is based entirely on the 
observed spatial symmetries of these patterns. However, information on spatial symme- 
tries of the new pattern alone may not be sufficient for our approach to be applicable 
in some problems. For example, consider a bifurcation problem defined on a spherical 
domain. Suppose a basic state with 0(3) symmetry loses stability and the observed 
bifurcating solutions are axisymmetric, then the isotropy subgroup of the bifurcating 
modes is given by 0(2). If the eigenfunctions are expanded in spherical harmonics, it 
is known that 0(2) is a maximal isotropy subgroup of 0(3) for all even values of the 
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spherical harmonic index I [14] and so an infinite number of irreps is relevant to the 
observed bifurcation. This example illustrates the fact that our method breaks down if 
the observed symmetries of the bifurcating modes form an isotropy subgroup for more 
than one irrep of the symmetry group of the basic state. 

We are currently involved in applying a similar method to the study of the 'superlattice- 
one' pattern reported in [1] as a bifurcating instability from standing hexagons. Pre- 
liminary analysis of the experimental data reveals that a suitable periodic box in this 
case will give rise to an arrangement of standing hexagons with a 'hidden' reflection 
symmetry [13], which leads to extra complications in determining the spatial symmetry 
group. It is an interesting problem that deserves further investigation. 

Unlike some time-periodic solutions (for example, standing rolls), which can also be 
defined on a hexagonal lattice, standing hexagons possess only trivial spatio-temporal 
symmetries [26]. So our treatment of the superlattice patterns as symmetry- breaking 
instabilities from standing hexagons is relatively simple because only instantaneous spa- 
tial symmetries are needed to define the isotropy subgroup of these solutions. In general, 
our approach can be applied to the study of spatial period-multiplying bifurcations from 
solutions with spatio-temporal symmetries and used to investigate some of the possible 
symmetry-breaking behaviour, if techniques discussed by Rucklidge & Silber [25] and 
Lamb & Melbourne [28] are also included. 
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